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Quasiclassical quantization of crystal dislocations field is considered in terms of func- 
tional integral. The generalized zeta-function is used to evaluate the functional integral 
and quantum corrections to mass in quasiclassical approximation. The quantum correc- 
r~| I tions to few classical solutions of one-dimensional Sin-Gordon model are evaluated with 

account of rest n — 1 dimensions. The results are applied to appropriate crystal dislocation 
. . models. 

1 Introduction 



^ ■ We consider Frenkel-Kontorova models for crystal structure dislocations [H [2] and its Sin- 
. Gordon (SG) equation counterpart. Being a one-dimensional field theory, which continuous 
limit is based on nonlinear Klein-Fock-Gordon (KFG) equations, it should be embedded into 
^ I two- or three-dimensional picture with effective account of rest variables as in the model itself 

as in quantization procedure. 
O ' The history of kink/soliton field quantization, started from [31 HI |5] is rather impressive 

now [6l [22] . One of advanced method, that allows to obtain relevant results in compact form 
is reviewed in [5], the semi-classical quantization is a well known method of approximating 
, energies of nonlinear equation's states. 
• Beginning from Maslov paper [7| the functional integral method becomes practical tool for 

evaluation of quasiclassical corrections to the action. At the meantime in [3] and in |1] such ex- 
pressions were studied for SG and (j)'^ models: approximate quantum corrections were evaluated. 
The problem of embedding such model into real multidimensional theory and the problem of 
regularization is still under examination [5] [ID], [12] • Generalization for the supersymmetric 
kink is given in 

Aim of this work is to apply the semi-classical quantisation method to calculate energy 
correction for Sin-Gordon kinks in real crystals. A general algebraic method of quantum cor- 
rections evaluation based on zeta-function [8] is used and the Green function for heat equation 
with soliton potential is constructed. The mentioned problems of the theory dimensions and 
regularization are traced. 
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In the first section connection between Frenkel-Kontorova model for a real solid and Sin- 
Gordon equation will be given and the kink solution of Sin-Gordon equation is reproduced with 
the realistic parametrization scheme. In the next section a space-time approach to the continual 
integral evaluation [23] and its semi-classical approximation of energy for static field equation's 
solutions will be described as well as the zeta-function regular isat ion method. Following section 
will be devoted to calculating corrections to energy of Sin-Gordon kink with expressions via 
constants of some solids. In last two sections these corrections are linked to Frenkel-Kontorova 
model of dislocations and crowdions. 



2 Frenkel-Kontorova model and Sin-Gordon equation 

Frenkel-Kontorova model describes one-dimensional chain of atoms interacting with each other 
and with a sinusoidal potential. Motion equation can be written as [2] 

g^x- dV 
M^ = G{xi+i-2xi + Xi_i)-—{xi), (2.1) 

where V (the external potential) usually takes the form 

\/(x) = |fl-cosf^')V (2.2) 



2 V V « . . 

where M is the mass of atoms of which the particular crystal lattice is built, a is the lattice 
constant, G is proportional to the bulk modulus K oi & given material (G = Ka) and e is 
proportional to shear modulus Mg {e = ^^^^ - only interaction with nearest neighbors is 
taken into account). It is also important, to note, that Xi can be viewed as actual position of 
particular atoms or deviation from equilibrium position (if equilibrium distance between atoms 
is equal to the external potential's period). The latter interpretation will be useful for obtaining 
continuous approximation (Sin-Gordon equation). Let w be a function of real variable x and 
let it solves the equation (12. ip on condition: 

u{ia,t) = Xi, Vz G Z. (2.3) 

By expanding all instances of function m in a Taylor series around ia and omitting higher 
derivatives of u, one obtains Sin-Gordon equation 

For the purpose of this work, we will rewrite the equation in unitless variables, with T as a 
parameter connected with action S (specified in the next paragraph): 

X = ax' , (2.5) 

u = a(f), (2.6) 

t = Tt'. (2.7) 

That yields 

|l||-<^0-5=in(..,)^O. (.8, 



Note, that the change of sign in front of cosine in the potential only moves the solution by 
a constant, thus doesn't have any notable effect. Well known static elliptic solution [12] of 
Sin-Gordon equation can be expressed as 



(p = — arcsin (A; sn {mx'; k)) H — . (2-9) 

TT 2 

with m = ^/^^TT, k G [0, 1]. Kink case is obtained for k = \. Classical kink energy in notation 
of 

= -VSea^G. (2.10) 

TT 

3 Method examination 

Let us consider a one- dimensional nonlinear Klein-Gordon- Fock (KFG) equation with t' as time, 
x' as spatial coordinate and a twice-differentiable potential ^ (1 — cos {2Tnp{x'))) 

It is the Euler equation of the variational principle with the action defined as 



S{ip) = Ta / dt' / dx' 

Jo J-oo 



aM [d^{x',t')V aG [ d^{x' ,t')V e 



r^^ i 1 1 r^ 1 [l — COS [211 LOix , t' ))) 

2r2 \ dt' J 2 \ dx' J 2a ^ v v-v . 

(3.12) 

Feynman formulation of quantum field theory links classic systems with their quantum coun- 
terparts through usage of the classical action in the path integral defining the propagator 



(V'le-i^^l^) = / D0(x',t')e-^^^*) (3.13) 



^0,1 



with H as quantum Hamiltonian and Cg'^ as family of all continuous functions fulfilling bound- 
ary conditions: 0(x',O) = 6{x'), 0(a;', 1) = tfj^x'). Let us consider a diagonal element of the 
propagator for a static solution ip 

{ip\e-^™\ip)= [ D0(x',Oe-^^(*^ (3.14) 



For further calculations, it is convenient, to represent as 

(p = if + '^ai(i)i (3.15) 

i 

with 0j as base vectors, that vanish at t' = and t' = 1 to satisfy boundary conditions. The 
Feynman integral can be evaluated by the use of stationary phase method [7] (since any classical 
solution is a stationary point of the action functional) 

= + + (3.16) 

j,k 



Second derivative of S can be expressed as 

^(^) = ra(^^„(^-— — + .G^-— co=(^2x(^ + X:..*.)jj^,j (3.17) 

with scalar product defined as 

(0fc,0,) = ^ dt' j dx'<Pu<P, (3.18) 

witli limits on the spatial integral to be specified depending on particular problem. Inserting 
it back to the propagator (13.141) gives 

(^le-^^^ly.) ^ e"*^(^) I D0(x',t')e^^^-^-'^^"'=^'^(^"("^^^'^''^"'^'°^^'"^^)'^^). (3.19) 

""-■0,0 

It is important to note, that static solutions of the field equation are assumed to be eigenstates 
of the Hamiltonian (or adequate eigenstates are approximated by pure classical states) 

e-i^^'(^|(^)^e"t^(^) ! D0(x',Oe-^^-^-"^"^(^^-'(-'^^^'^'''^^-'^™^^'^^^)'^^), (3.20) 

'-'0,0 

where Eq is the quantum energy of the state Lp. Since —^^-^ + a^GT-^ — 2£7r^Tcos (27r(/j) 
is a sum of one-dimensional second-order linear Hermitian operators, it's eigenf unctions form 
a complete set of orthogonal functions. In this situation one can choose the set {0j} to be 
composed of said eigenf unctions 

e-'nTE.{^\^) ^ e-i^(-) / D0(x',t')e-^^^'^K<^-(-^^+'^«^^-^^-M<^O. (3.21) 

7^0,1 



Also D(j){x' ,t') can be rewritten as [2Tj 

D(t){x' ,t') = \[dak, (3.22) 

k 

nrfa.e-^^^'^K^-(-^^+'^«^^-^™^M^O. (3.23) 

■°° k 

From orthogonality one obtains 



nrfa.e-^'^'^('^-(-^^+'^^^^-'^™^M^O. (3.24) 
k 

If Aj will be considered as eigenvalues of — + ^'^^T-^ — 2£7r^rcos (27r(/?), then 

/oo 

JJrfafcC-^'^'-^^^'^^''^^). (3.25) 
-oo 



Substitution: 



A 



k {(Pk,(Pk) 

2h 



2h 



Xk 



k {(Pk, Vk) 



2h 



Xk 



k \ (pk, (pk 



da',,e-<. 



The remaining integral gives: 



J — c 



da'^e-'^' 



2nh 



iXk {(t>k, (t>k) 



As for the norm of (pk 



5fc, <Pk 



dt' I dx'4)\, 



(3.26) 
(3.27) 

(3.28) 

(3.29) 
(3.30) 

(3.31) 



it will be considered unitless with norm of 1, which is ensured by previous change of variables 
^2M . Thus: 



e ''{(f\ip) ~ e " 



'n 



2Txh 
iXk ' 



e-^^^'((^|(^) ~ e-*^(^ydet 
Eg can be now expressed as 



27Th 



+ a^GT^ -2en^T cos {2n^) 



For brevity 



2iT 



In det 



27ih 



(3.32) 



(3.33) 



D 



2Txh 



-^W2+-'GT£,^-2e.^Tco.{2.^{x')) 



(3.34) 



(3.35) 



Problem: differential operators have energy spectrum up to infinity, thus their determinants 
will be infinite as well. First step in solving the problem, is to properly set zero for the energy 
by subtracting an analogous term for vacuum [12] 



h 



h 



h 



Eq ~ E{ip) + — In det [D] - — In det [Do] - -7^\n{ip\(p) 



2iT 



with 



2iT 



iT 



(3.36) 



Second step is to rewrite the determinants in a form that allows subtraction. For this purpose, 
let us define generalised zeta-function as 

c,,(.) = Yl (3-38) 

fc:Afc>0 

where are eigenvalues of a given operator D. One should note, that summation is done over 
eigenfunctions, also the summation itself should be treated in the most general sense. It is 
evident, that [12] 



dC 

lndet[D] = -^(0). (3.39) 

To obtain C functions, when spectra of involved operators are unknown, one should consider 
an equation 

■^ + D^ go {y, x', x'o) = S{y)d {x' - x'q) . (3.40) 
Green function for this equation can be written as (for y > 0) 

gn {y, x', x',) = J2 {x') rk {x'o) O (y) , (3.41) 

k 

where ipk are orthonormal eigenfunctions of operator D. 

Only positive eigenvalues of the operator D are important for further calculations, so one 
can integrate over all spatial variables (the resulting function is defined as 7 function) and 
subtract e°* (0 is an eigenvalue for any linear operator): 

Id (y) = -1 + J go (y, x, x) dx, 

ID (y) = Yl (3-42) 

fc:Afc>0 

By performing Mellin transform 

r+00 

Cd{s) = tttt y"'^iD{y)dy, 
r(s) Jo 

one obtains generalized Riemann zeta-function 

Us) = (3.43) 

fc:Afe>0 

It is important to note, that renormalization done in f l3.36p is not complete. There as 
yet arises an ultraviolet divergent term [3], which can be countered by multiplication of the 
zeta-function by an appropriate term [12] 

„2s 



Cr{s) = r'X{s), (3.44) 



which, from a point of view of the definition fl3.43p . may be interpreted as the spectrum rescaling 
Afc r~^Afc. Thus 

2^(C'(0) + 21n(r)C(0))-A 



^ E{^) - — (C'(0) + 21n(r)C(0)) - — ln(v9|^). (3.45) 



For an operator that can be written as a sum of operators dependant on different variables 
D = Di + D2, J function is a product of 7 functions of the summed operators [12] 



7D = 7D^lD2, (3.46) 
what greatly simplifies calculations. In this work: 

Di = ^ (a^GT^ - 2eiT'T cos (27ryp(x'))) , (3.47) 

D2 = 7rT{ • (3-48) 



2nh V T dt'^^ 

4 Corrections for a Sin-Gordon kink 

We can now proceed to calculating energy corrections. For brevity let's define 

U{x') = -— cos {2TTip{x')) (4.49) 



as well as: 



t/(x') = 5^ (2*' - 1 - cn' ( J^^^; k I I (4.50) 



2'KhT' ^ ' 



m = \l^n. (4.53) 
Having the potential U, one can find the Green function for the operator Di. 

— + A— - AU {x)\ {go, {y, X, xo)) = 6{y)S {x - Xq) . (4.54) 

After rescaling yA = Ay we have 

+ ^ - {^qd^ (^^,x,xo)) = 5{yA)5 {x - xq) . (4.55) 

To obtain the Green function, one uses Laplace transform 

p + — - [/(x) ]gD^ {p, X, Xo) = 6{x-Xo), (4.56) 



dx 



where 



9di (J^,x,Xo^ = J^gn^{p,x,Xo)ePy^dp. (4.57) 



For further calculations one needs only values of for In such a situation 

gDi{p,x,x) = G{p,x), where G{p,x) solves equation, similar to Hermit one [T5].[T2] 



2GG" - {G'f - 4(t/(x) - p)G^ + 1 = 0. (4.58) 
To solve equation (14.581) . one substitutes 

z = cn^{mx] k), (4.59) 

which gives: 



2G{G'^,4m^{-ez^ + (2^ - l)z^ + (1 - e)z) + G',2m^{-3ez^ + {Ak^ -2)z + l- k^))- 
Am''{-k^z^ + (2A;2 _ ^ _ p)^) ^Qr^2 _ ^^^^^^ _ ^^^2 + 1 = 0, 



One postulates solution in a form 



(4.60) 



G{p,z) = ^^^^^, (4.61) 

where P i Q are polynomials: 

P = p + az + b, (4.62) 

Q = QsP^ + q2p'^ + qip + qo- (4.63) 

After separating the equation with respect to monomials of p, one obtains: 

P{p, z) = p - m^k'^z, (4.64) 

For the kink case {k = 1) function G{p,x) takes the form 

Gn, {p, x) = - ^^-66) 

2^171^— p 2py'm'^~p 

Solution for a constant potential (vacuum) can be obtained from f l4.60p with constant u{z) 

Gi,,{p,x) = -jl^. (4.67) 



At this point one can perform subtraction from ([X 



G..-G.„^G.(p..)^- f7";''"^' . (4.68) 

2p^ — p 



Integrating with respect to x one obtains 



Inverse Laplace transform of the above function can be found in mathematical handbooks: 
7i [-j^ ) = Erf{imy/yA) = J e ^ dx = — -^j^ J exp[m yAT Jdr, (4.71) 

71(2/) = f exp[m^Ayr^]dT. (4.72) 

Now one has to incorporate 'JD2 calculated below: 



7(?/) = 7i(y)7z)2(2/)> 



l{y) = liivhoiiy) = """"^ I exp[m^AyT^]dT. (4.73) 



l{y) = iMidAv) = 7^ / eicg[-nn?{-A)yT^]dT. (4.74) 

TTV B 



After applying Mellin transform: 



TiT{s) V S 7o 



+00 /•! 



C(^) = -^a/^ / 2/^"^ / e-^(-^)^^^dr, (4.75) 







Above integral converges if Re{s) < \. With such an assumption one obtains 



m^^'+H-A)-' -A 1 



Derivative of Q function: 



C'(0) = 7 V ^ (2 - 2 ln(m) - H-A)) . (4.77) 



Inserting into fl3.45p gives 



E, ^ a^G^ - y (2 - 2 ln(m) - In(-A) + 2 ln(r)) - A (4.78) 

After substituting m, ^r, A and B 

(4.79) 



The choice of the scahng parameter r allows to eliminate a dependence of the result on T (the 
simplest choice is = ^). 



This result for the correction coincides with one of [3]. 

5 Dislocations 

5.1 General scheme 

An outline for modeling edge dislocations through Frenkel-Kontorova model consists of two 
steps [2]. Firstly one simulates cross-section of the dislocation (in direction of Burgers vector) 
as a Frenkel-Kontorova kink. Parameters for the equation 



are obtained as follows: M is the mass of atoms of which the particular crystal lattice is built, a 
is the lattice constant, G is proportional to the bulk modulus K of a given material {G = Ka) 
and e is proportional to shear modulus Mg {e = ""^^ ) • 

At this point one can calculate energy of the dislocation and determine the Peierls-Nabarro 
potential [H], which is accountable for interactions between moving kinks and crystal lattice. It 
is important, not to approximate Frenkel-Kontorova equation with Sin-Gordon equation at this 
point, or else Peierls-Nabarro potential will vanish. Approximate solutions of Frenkel-Kontorova 
equation are often constructed on basis of Sin-Gordon soliton though. 

In the second step one describes the dislocation line as a row of Frenkel-Kontorova kinks that 
interact harmonically with their nearest neighbors and are subject to Peierls-Nabarro potential 
(or it's first term in Fourier expansion). This means, that the dislocation line is described by 
another Frenkel-Kontorova equation, which can be well substituted by Sin-Gordon equation. It 
is important to note, that the harmonic interaction coefficient as well as mass is not the same 



Movement of a dislocation is simulated through creation and propagation of second order 
kink-antikink pairs as a response to applied stress. Often the above described framework is 
enriched by taking into account additional terms of the potential or anharmonic interactions. 
It might also be important to consider thermal oscillations of the lattice in modeling macroscopic 
properties of crystals. 

5.2 Calculations 

Peierls-Nabarro potential's amplitude (^2) and harmonic interaction coefficient {G2) for 
dislocation line modeling were calculated numerically. Firstly approximate static Frenkel- 
Kontorova kinks were obtained by substituting Sin-Gordon kink as an initial condition to a 




(4.81) 




(5.82) 



as in fl5:82|) . 



critically damped Frenkel-Kontorova equation (already in finite difference form) 



fixi, tj) = f{xi, tj^i) + {f{xi+i,tj^i) - 2f{xi, + f{xi^i,tj^i))dt+-^ sm{27i f{xi, tj^i))dt. 

(5.83) 

By changing initial position of the kink, one gets two different stable Frenkel-Kontorowa 
kinks, which correspond to minimal and maximal kink energy. Difference of their energy gives 
£2- G2 is estimated on assumption, that energy of any two atoms interacting is given by 

E{r) = ^ - y > (5-84) 
where Ci and C2 are fitted so that E(a) is the minimum and E"{a) = G 

Gi = \a^G , . 

G2 = a'G ■ ^^-^^^ 

Interaction energy for a given displacement between kinks (AX) is then calculated by adding 
up interaction energy of atom pairs from two neighbouring kinks. To avoid major rounding up 
inaccuracy and set in the zero for energy, one subtracts an analogous term for AX = on every 
iteration. By varying AX, one can obtain G2 by fitting ^AX^ as the energy. It is important, 
to fit this function separately for both Frenkel-Kontorova kinks, which means that AX = na. 
Effective mass of the kink is obtained thanks to the fact, that Sin-Gordon equation is Lorentz 



invariant (with sound velocity c = J used instead of speed of light in vacuum) 



M. = ^ = zmJ§^. (6.86) 



where Eq is the energy of a static kink. 



5.3 Application 

Aim for this subsection is to estimate energy corrections to a kink-antikink pair approx- 
imating it with a sum of kink and antikink (which is acceptable when kink and antikink are 
far enough from each other). This in turn will have effect on calculated plasticity of a given 
material [16]. As pointed in section (15. ip all important coefficients can be obtained from bulk 
modulus, shear modulus and lattice parameters for a given crystal. It is to be noted, that one 
can only estimate energy correctly for pure elements monocrystals. Further calculations in this 
subsection are based on parameters of polycrystal metals (due to lack of more relevant sources), 
thus are only estimation of order of quantum corrections for kink energy in dislocations. As 
can be seen in Tabled] [E^ as classical energy and AE^ as corrections), corrections are three 
orders of magnitude lower then the classical energy. Therefore they should be considered if 
desired precision is of three significant figures or more. 



6 Crowdions 



The most basic model of a crowdion [2] is similar to the first step of dislocation modelling. 
Essentially the same Frenkel-Kontorova equation is used with the exception, that the potential 



now describes interaction with all atoms surrounding the modelled chain of atoms. For the 
purpose of this section the effect of external stress on crowdion energy will not be included as 
well as deformation of the lattice around the crowdion (on contrary to [19] or [20] ) 



M 



G (xj_|_i — 2xi + X 



i-l) 



£7r . .27rxj, 
— sm( 

a a 



(6.87) 



where M is the mass of atoms of which the particular crystal lattice is built, a is the, G is 
proportional to the bulk modulus of a given material (G = Ka) and e is proportional to 
shear modulus Ms {e = ^^^^ - only interaction with nearest neighbors is taken into account). 
The above equation is approximated by Sin-Gordon equation as described in section (|2|) 



Md'^u 



aG 



dx"^ 



en . 
+ ^sm 



2nu 



0. 



(6.88) 



Crowdion classical energy and quantum corrections are calculated as shown in sections (|2|) and 





M [kg -10-^"] [S] 


a [nm] |17| 


Ms [GPa] [18] 


K [GPa] [18] 


Ed [meV] 


AEa [meV] 


Ec [eV] 


AEc [eV] 


Ag 


17,9119 


0,40776 


30 


100 


1,9229 


0,0113 


9, 3934 


0, 0034 


Al 


4, 48345 


0,40494 


26 


76 


2, 7443 


0, 0336 


7, 4665 


0, 0064 


Au 


32,7071 


0,40702 


27 


180 


0,06719 


0, 00033 


11,8909 


0, 0024 


Cu 


10,552 


0,3610 


48 


140 


3,1938 


0, 0283 


8, 7556 


0, 0053 


Fe 


9, 2733 


0,28665 


82 


170 


7, 4299 


0, 0996 


7, 0290 


0, 0066 


Mg 


4,03594 


0,5196 


17 


45 


5, 5603 


0, 0449 


9, 8149 


0, 0061 


Ni 


9, 7463 


0,429 


76 


180 


14,4136 


0, 0786 


23,3596 


0, 0076 



Table 1: Quantum corrections to dislocation kink and crowdion energy for exemplary metals 



Corrections (see Table Hj) are four orders of magnitude lower then the classical energy, yet 
they can be notable in precise measurements. 



7 Conclusion 

An important element of the approach which we demonstrate in the Sec. 4 relates to the 
renormalization problem. Account of the "vacuum" contribution with the constant potential 
term gives the finite result for the energy correction (compare with e.g. [22]). 
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